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Abstract 

We construct a Gaussian approximation to the effective theory for the Colour 
Glass Condensate which describes correctly the gluon distribution both in the low 
density regime at high transverse momenta (above the saturation scale Qs), and 
in the high density regime below Qg, and provides a simple interpolation between 
these two regimes. At high momenta, the effective theory reproduces the BFKL 
dynamics, while at low momenta, it exhibits gluon saturation and, related to it, 
colour neutrality over the short distance scale 1/Qs ^ 1/Aqcd- Gauge-invariant 
quantities computed within this approximation are automatically infrared finite. 



1 Introduction 



There has been significant progress towards the understanding of the high energy, or 
small-x, limit of QCD in the recent years, and part of this progress is associated with 
the construction of an effective theory for the small-x gluons in the hadron hght-cone 
wavefunction: the effective theory for the Colour Glass Condensate (CGC) |], ^ |^, |^, 
1^, H, H (see also Refs. |^ |ll| for recent reviews and more references). This is an effective 
theory since it applies to gluons with a given value of the longitudinal momentum fraction 
X, with X -C 1, and is obtained after integrating out the gluons with larger values of x 
within QCD perturbation theory. 

The general strategy behind this construction is reminiscent of the Born-Oppenheimer 
approximation (i.e., there is a separation of scales between "fast" and "slow" degrees 
of freedom), and the resulting mathematical description is that of a glass. The "fast" 
degrees of freedom are the gluons (or, generally, partons) with larger values of x, which 
do not participate directly in the high-energy scattering (because of their large relative 
rapidities), but act as sources for the small-x gluons which produce the scattering. The 
"slow" degrees of freedom are the small-x gluons themselves, which have low rapidities 
(of the same order as the external projectile), but large densities in the impact parameter 
space (the transverse plane), since their radiative production is enhanced at small-x. 

This whole partonic system is characterized by weak coupling — since the high gluon 
density acts as a hard scale for the running of the QCD coupling — , but strong interac- 
tions with the small-x gluons, because these latter have densities which are parametrically 
of order l/ctg. The crucial simplification, which allows for the all-order resummation of 
the density-enhanced interactions, is that the small-x gluons can be treated in the clas- 
sical approximation [Q, because of their large occupation numbers. That is, the small-x 
gluons can be described as the classical colour fields radiated by fast moving "colour 
sources" which represent the partons with larger values of x. Then, the all-order resum- 
mations alluded to before amount to treating the non-linear effects associated with these 
classical fields exactly , both in the classical field equations [|I|, 0, |^, Q , and in the quantum 
calculations involved in the construction of the effective theory |p, |^. 

A second simplification, coming from the separation of scales between "fast" and 
"slow" , is that the internal dynamics of the fast partons can be neglected over the char- 
acteristic time scales for the dynamics at small x. This is the premise of the "glassy" 
description : For the purposes of the small-x physics, the fast partons can be treated as a 
frozen configuration of classical colour sources, which is random, in the sense that it can 
be any of the configurations allowed by the dynamics of the fast partons, and which must 
be averaged over in the calculation of physical quantities. This averaging is necessary also 
to restore gauge symmetry: on the average, the colour charge density must vanish at any 
point. This means that the weight function for this averaging, i.e., the probability density 
for having a given configuration of the colour sources, must be a gauge-invariant func- 
tional of the colour charge density. Further characterizing this functional, and deriving 
an explicit expression for it, will be our general objective in this paper. 
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This weight function encodes all the information about the colour source and its cor- 
relations, as inherited from the quantum dynamics at larger values of x. When x is further 
decreased, the structure of the classical field equations does not change: these are always 
the Yang-Mills equations with a colour source in their right hand side. What do change, 
are the correlations of the colour source, that is, its functional weight function. This 
change can be computed in perturbation theory, and expressed as a functional renormal- 
ization group equation (RGE) for the weight function P, 0. The initial condition for this 
equation is non-perturbative, and requires a model. But once this is specified, then the 
effective theory at scale x is completely determined by the solution to the RGE. 

Even though the general solution to this equation is not known (this is a complicated 
non-linear equation 0; see, however, |jl2[), its physical implications are well understood, 
at least, in some limiting kinematical regimes, where approximate solutions have been 
obtained P, [1^. One has thus identified two important physical regimes: (I) a low- 
density regime at high transverse momenta (for fixed x), where the non-linear effects can 
be expanded out, and the RGE reduces to the BFKL evolution |jl3[, and (II) a high- 
density regime at low transverse momenta, where the non-linearities are essential, and 
the RGE predicts gluon saturation p!5| , |16|. The transition between these two regimes 
occurs at transverse momenta of the order of the "saturation scale" Qs{x), which is the 
typical momentum of the saturated gluons, and grows like a power of 1/x Jl^, [1^, |19|, . 

and numerical [|^, 



These conclusions are corroborated by analytic pi] , P2| , |T8 



23], studies of the Balitsky-Kovchegov (BK) equation which describes deep 

inelastic scattering at high energy, and is encoded too in the RGE 0. Originally, the 
BK equation has been derived in a different approach |T^, ^ (see 

also Refs. ^ for earlier versions of this equation), which focuses on the observables 



through which an external projectile couples to a hadronic target, rather than on the 
wavefunction of the target itself. 

Our purpose in this paper is to use the conclusions of these previous studies to con- 
struct an approximate solution to the RGE which is simply a Gaussian in the colour 
charge density p. Clearly, having an explicit expression for the weight function, which 
is moreover a Gaussian, will greatly facilitate the use of the CGC effective theory in 
practical calculations. In fact, with such a weight function, most calculations within the 
effective theory can be done quasi-analyticallyQ. Besides, a Gaussian weight function can 
be easily implemented in numerical simulations, like those performed in Ref. ||3^, where 
the effective theory is used to study nucleus-nucleus collisions on a lattice. 

The adequacy of a Gaussian approximation may look at first surprising, given the 
highly non-linear structure of the RGE. By definition, a Gaussian approximation retains 
only the 2-point correlation function (pp) among all the ra-point functions of the colour 
sources. Thus, this is a reasonable approximation provided (a) one disposes of a good 
approximation for the 2-point function, and (6) the higher n-point correlations are rela- 



^Note that, even with a Gaussian weight function, the CGC effective theory remains non-trivial, since 
the classical field equations (i.e., the Yang-Mills equations with colour source p) are non-linear. Still, as 
we shall remind in Sect. 2, the exact solution to these equations is known in analytic form. 
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tively weak. But in the RGE, the evolution of the 2-point function is generally mixed with 
that of the higher n-point functions, so it is not a priori clear how to construct a sensible 
approximation for the 2-point function alone (and thus a Gaussian weight function). It is 
even less clear whether such an approximation can be relevant for the non-linear regime, 
where higher correlations may be important as well. 

Still, as we shall demonstrate below, a Gaussian weight function can accommodate 
both the BFKL evolution of the gluon distribution at high transverse momenta, and the 
phenomenon of gluon saturation at low transverse momenta. It is only the transition 
between these two regimes that cannot be accurately described in this way. (But our 
approximation will provide a smooth interpolation between these limiting regimes.) To 
understand when, and to which extent, is a Gaussian approximation expected to work, 
let us make a few remarks that will get substantiated in the subsequent analysis: 

In the low-density regime at high transverse momenta k± ^ Qs{x), the evolution 
generated by the RGE does not couple correlation functions with different number of 
external legs. In particular, the 2-point function (pp) satisfies a closed equation, which 
is moreover linear @]. This is the BFKL equation |13|. By solving this equation, one can 
construct a Gaussian weight function which incorporates correctly the BFKL dynamics 
for both the unintegrated gluon distribution, and deep inelastic scattering. 

In the saturation regime at low momenta k± -C Qs{x), the colour sources and the 
radiated fields are typically strong (their amplitudes are parametrically of order l/g), so 
the dynamics is fully non-linear. Remarkably, however, the non-linear effects drop out 
completely from the RGE. This is related to the specific way how these effects enter the 
RGE in the first place [0] : they enter only via Wilson lines, that is, path-ordered expo- 
nentials of the classical fields. At saturation, where the fields are strong, these complex 
exponentials oscillate rapidly and average to zero, together with all their correlations, 
when probed over transverse distances of the order of the saturation length 1/Qs{x), or 
larger. Thus, in this high-density regime, the RGE becomes effectively quadratic^ and 
can be trivially solved. The resulting weight function is then truly a Gaussian 

We see an interesting duality emerging at saturation: this strong field regime allows 
for a description in terms of a Gaussian weight function, so like a free theory. 

Another important feature at saturation is colour neutrality p, |^, |10|, ^ : The satu- 
rated gluons shield each other in such a way that their global colour charge vanishes when 
integrated over a disk of radius 1/Qs{x), or larger. This is expressed by the following 
property of the charge-charge correlator: {p{kj_)p{—k±)) oc k1_ for momenta kj_ <C Qs{x). 
This property has crucial implications: 

First of all, it is synonymous of saturation: As we shall see, {p{k±)p{—k±))/k1_ is 
essentially the unintegrated gluon distribution. Thus, the fact that the charge-charge 
correlator is proportional to k'^ at low momenta, means that the unintegrated gluon 
distribution is independent of k^ (up to logarithms), which is saturation indeed. In fact, 

^In the same sense, for instance, as the Hamiltonian for a quantum-mechanical harmonic oscihator is 
a quadratic operator. 
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we shall see that the scale Qs{x) is primarily generated (by the non-linear effects in the 
quantum evolution) as the typical scale for colour neutrality. Then this becomes also 
the "saturation scale" because of the connection between the distribution of the colour 
sources and that of the small-x gluons alluded to above. 

Second, colour neutrality ensures that gauge-invariant quantities come out infrared 
finite when computed in the effective theory. This should be contrasted with the original 
McLerran-Venugopalan (MV) model in which the colour sources are assumed to be 
uncorrelated, so the corresponding charge-charge correlator {p{k±)p{—kj_)) is independent 
of k±. Because of that, the calculation of physical quantities like the gluon distribution 
function Q], or the dipole-hadron scattering amplitude |T^, within this model meets with 
logarithmic infrared divergences. It has been proposed by Lam and Mahlon to remove 
these divergences by imposing colour neutrality on the size of the nucleon, that is, by 
effectively taking into account the non-perturbative mechanism of confinement. In Ref. 
P^ , this suggestion has been implemented in numerical simulations, but the results turn 
out to be quite sensitive to the poorly known non-perturbative physics. As we shall show 
in this paper, when the colour neutrality due to saturation is taken into account, all such 
infrared divergences are eliminated already at the perturbative scale 1/Qs{x). 

Other effects which can be attributed to colour neutrality at saturation, and which 
illustrate the lack of infrared sensitivity of the effective theory, are the suppression of the 
"infrared diffusion" in the quantum evolution towards small x (this has been verified in the 
numerical simulations of the BK equation ) , and the predominance of the short-range 
scattering in the evolution of the total cross-sections with increasing energy 

Note that colour neutrality is a weaker property than confinement: Although the 
global colour charge vanishes indeed, the multipolar moments (starting with the dipolar 
one) are still non-zero, and generate long range interactions. But the latter die away suf- 
ficiently fast to ensure the infrared finiteness of all observables for high-energy scattering 
(at fixed impact parameter). The information about confinement is needed only in the 
calculation of the total cross-section, since this requires an integration over all impact 



parameters 0, Still, as shown in Ref. |^ on the example of the dipole-hadron 

scattering, by combining the non-linear evolution equations with a minimal assumption 
about confinement in the initial conditions, one can compute the rate for the increase 
of the total cross-section with the energy, and conclude that citot saturates the Froissart 
bound 1 3^ at very high energy. Similar conclusions have been reached in Refs. |3^, 37 . 



In view of the previous considerations, our general strategy in this paper should be 
clear by now: We shall construct our Gaussian approximation for the weight function 
by requiring its kernel (pp) to satisfy the BFKL equation at high momenta, to describe 
colour neutrality at low momenta, and to interpolate smoothly between these two limiting 
regimes, with the change in behaviour occurring at the saturation scale. The last condition 
requires a saturation criterion for the charge-charge correlator, that we shall formulate, 
and prove to be equivalent to the other criteria used in the literature. For practical 
purposes, we shall propose also a simplified form of the kernel, in which the general 
solution to the BFKL equation is replaced by its "geometric scaling" approximation , 
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valid up to momenta k±_ ~ Ql{x) / Kqcd- 

To test the quality of our approximation, we shall compute the unintegrated gluon 
distribution and the dipole-hadron scattering amplitude, and find the expected results 
at both high and low momenta. In particular, from the kernel of the Gaussian, we 
shall deduce a simple analytic formula for the unintegrated gluon distribution, which 
has a scaling form, i.e., it is a function of k_i_/Qs{x), and interpolates smoothly between 
saturation at k± < Qs{x) and BFKL behaviour at > Qs{x). 

Also, we shall derive a phenomenologically useful formula for the dipole-hadron scat- 



tering amplitude, which properly incorporates the constraints of geometric scaling 



(both in the saturation regime, and in the scaling window |T8[ above Qs{x)), and the 



colour neutrality due to saturation. This analytic formula can be easily applied to the 
description of electron-hadron deep inelastic scattering, as an alternative to the more 



phenomenological parametrizations proposed in Refs. ||39| , ^ |. 

At this point, let us note an important difference in the mechanism of saturation 
between the Gaussian effective theory that we shall construct and the original MV model 

which has a Gaussian weight function too: In the MV model, the colour sources 
are uncorrelated, and the saturation arises exclusively via the non-linear effects in the 
Yang-Mills dynamics of the classical fields. By contrast, in the effective theory which 
includes quantum evolution, the saturation is the consequence of the long-range spatial 
correlations between the colour sources, and, as such, it is encoded already in the weight 
function. This explains why, in our case, the unintegrated gluon distribution can be simply 
proportional to the kernel of the Gaussian even in the non-linear regime at saturation. 

This paper is organized as follows: In Section 2, we shall succinctly review the effective 
theory for the CGC, and introduce the notations and conventions to be used through- 
out. In Section 3, we show that Gaussian approximations for the weight function can 
be formally generated through mean field approximations to the RGE. In Section 4, we 
briefly describe the MV model ^ |^, ^, which will serve as a basis of comparison for 
the following developments. In Section 5, we discuss the simpliflcations which occur in 
the RGE in the limiting situations where the transverse momenta are either very large, 
or very small, compared to the saturation scale. This justifles the use of a Gaussian ap- 
proximation in these regimes, at least. In Section 6, we construct an approximation for 
the charge-charge correlator which interpolates smoothly between the limiting behaviours 
found previously. This specifles the Gaussian effective theory completely. In Section 7, 
we use this effective theory to flrst derive a simple expression for the unintegrated gluon 
distribution, and then compute the dipole-hadron scattering amplitude. Our conclusions 
are summarized in Section 8. 
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2 The Colour Glass Condensate 



2.1 The classical effective theory 

Let us summarize here the equations at the basis of the effective theory for the CGC. 
The classical field obeys the Yang-Mills equations with a colour current due to the fast 
parton^ : 

= 5^+pa{x-,x^). (2.1) 

The structure of the colour current in the r.h.s. is fixed by the kinematics: i) The 
current has just a plus component (J^ ~ S^^pa) since the fast partons move nearly at 
the speed of light in the positive z (or positive a;"*") direction, ii) By Lorentz contraction, 
the support of the colour charge density pa is concentrated near the light-cone, i.e., near 
x~ = 0. in) By Lorentz time dilation, pa is independent of the light-cone time x~^, i.e., 
it is "frozen" , according to the terminology used in the Introduction. 



All the approximations involved in writing eq. (2J.) are consistent with the quantum 
evolution towards small-x in the leading-log approximation, since this evolution preserves 
the separation of scales in longitudinal momenta (and therefore also in time). 

The observables that can be computed in the effective theory are the gauge-invariant 
correlations of the gauge fields. These are first evaluated on the solution A^[p] to eq. (p.l|), 
and the result is then averaged over all the configurations of p, with a gauge-invariant 
weight function W^t[p] which depends upon x via the rapidity variable r = ln(l/x). 

For instance, we shall need later the scattering amplitude for the high energy scattering 
between a "colour dipole" (say, a quark-ant iquark pair in a colourless state) and the 
hadron. This quantity is a physical observable (at least, indirectly) in the sense that 
it enters the formula for the F2 structure function at small x p7|, |4l|. In the eikonal 



approximation, the scattering amplitude is obtained as ||42| , |26| : 

J^T{x±,y_L) = 1 - Sr{x±,y±), Sr{x±,yi_) = (tT(v\x_L)V{y_L)^^ ^, 

(tr(rt(x^)V(y^)))^ = Jdp WAp] tT(vl[p]VyAp]), (2.2) 

where V^{x±) and V{y±) are Wilson lines describing the interactions between the fast 
moving quark, or antiquark, from the dipole and the colour field in the hadron: 

V\x±) = Pexpji^ j rfx" A+(x",xx)t"| , (2.3) 



•^We recall that, in light-cone coordinates, and with our present conventions, = {t + z)/^/2 plays 
the role of the light-cone time, while x~ = {t — z)/^/2 is the longitudinal coordinate. Correspondingly, 

is the light-cone energy, and is the longitudinal momentum; thus, x = k^ /P^, where fc+ refers to 
a parton in the hadron, while denotes the total longitudinal momentum of the hadron. 
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and = A+[p] is the projection of the classical colour field along the trajectory of 
the quark (or the antiquark). The Wilson line correlator in eq. ( |2.2| ) is written in a 
suitable gauge in which the component is non-zero, but the result is actually gauge- 
independent, since this correlator can be completed to a Wilson loop by joining the points 
x± and y± via Wilson lines in the transverse planes at x~ = ±00 (where the transverse 
fields y4* can be chosen to vanish in any gauge with 7^ 0). 

Explicit calculations are most conveniently performed in the covariant gauge 9^^^ = 
0, since in this gauge the solution A'^ to the classical equation ( p.l|) has only one non-trivial 
component: A'^ = 6'^~^aa, with satisfying the two-dimensional Poisson equation: 

-V]_aa{x~,x±) = Pa{x~,x±). (2.4) 

This simplification reflects the fact that the classical field radiated by p has only one 
non-trivial field strength, namely the "electric field" which in the covariant gauge 
is computed as —d^aa- The corresponding result in other gauges is then obtained via the 
appropriate gauge rotation. Particularly interesting for what follows is the result in the 
light-cone (LC) gauge ^4+ = 0, which reads (in matrix notations; see J^, |10| for details) : 



F+\x-,x±) = U{- d'a) , (2.5) 

where 

U\x~,Xi_) = Pexp-cigf / dz~ a{z~ , X i_) > , (2.6) 



is the matrix for the gauge rotation from the covariant gauge to the LC gauge. This is 
essentially the same Wilson line as for the quark scattering in the eikonal approximation, 
eq. (p73|). Note that eqs. (|2.4|) - (|2.6| ) provide an explicit expression for the electric field in 
the LC gauge in terms of the colour source pa in the covariant gauge. In this sense, the 
solution to the classical field equations ( |2.1| ) is known exactly (in any gauge). 
Eq. (p.4|) can be easily inverted : 



aa{x ,Xi_) = I d'^y±{x±\—:^\y±) paix ,y±) 

The infrared cutoff A is necessary to invert the Laplacian in two dimensions, but, as we 
shall see, it disappears in the calculation of physical quantities. From eq. (|2.7|) , it is clear 
that the "Coulomb field" a{x~,x±) is time-independent and localized near x~ = 0, like 
the colour source. 

To summarize, to compute the dipole-hadron scattering amplitude ( |2.2| ) within the 
effective theory, one has to replace ^4+ aa in eq. ( |2.3| ), with given by eq. (|2.7| ), and 



then perform the functional integral over p in the second line of eq. ( |2.2| ), with the weight 
function Wr[p] (which for this purpose is needed as a functional of the colour source pa 
in the covariant gauge). 
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2.2 The Renormalization Group Equation at Small x 



When decreasing x (or increasing the rapidity r = ln(l/x)), new quantum modes become 
relatively "fast", and must be included in the effective source at the new value of x. 
This operation modifies both the longitudinal support of p (with increasing r, the colour 
source extends to larger and larger values of the x~), and its correlations. All these 
changes can be absorbed into a functional renormalization of the weight function PVV[p]; 
which therefore depends upon r. To the order of interest, this dependence is governed by 
a functional renormalization group equation, which reads : 

with x[p] a positive definite functional of p, whose explicit form is not needed here (this 
can be found in Refs. [0, |, ^^l)' but whose relevant properties will be discussed below. 

Eq. ( p.8|) allows one to construct the weight function by integrating out the quantum 
gluons in layers of r. The kernel x[p] encodes the changes in the correlations of the colour 
source induced at one step in this construction, and depends upon the source p created in 
all the previous steps. Indeed, the quantum gluons which are integrated out in one step 
propagate through the background field A[p] of the colour source generated in the previous 
steps. In general, this classical field is strong, so its non-linear effects on the dynamics 
of the quantum gluons must be included exactly. Because of that, the functional x[p] is 
non-linear in p to all orders, and also non-local, both in the longitudinal coordinates, 
and in the transverse ones. In fact, the non-linearity in p and the non-locality in x~ are 
strongly correlated, since they have a common origin: namely, the fact that x[p] depends 
upon p via the Wilson hues in eq. ( p.3|) . 

There is one more notation in eq. ( |2.8| ) which waits for an explanation: this is the 
argument p'^{x±) of the functional derivatives there. This is a notation for the colour 
charge density p"'{x~,x±) at a specific longitudinal location x~ = x~, which is the upper 
boundary of the support of the classical source in x~. That is, the colour source at rapidity 
T has support within a limited interval in x~ , namely at < < with x~ oc e^. 

This correlation between the quantum evolution in r and the longitudinal extent of the 
colour charge distribution can be understood as a consequence of the uncertainty principle: 
Since obtained by integrating out quantum modes with large longitudinal momenta ^ 
xP^ = e~^P^, the classical source at rapidity r must be localized near x~ = 0, within a 
distance Ax~ ~ e'^x^ (with Xq = 1/P+). 

But eq. ( p.8|) shows that this correlation is even stronger: Namely, it shows that, when 
the rapidity is further increased, say from r to r + dr, the additional contribution to the 
colour source which is generated in this way has no overlap in x~ with the original source 
at rapidity r. Rather, this new contribution makes a new layer in x~ , which is located 
between x~ and x~_^^^. This is why the functional derivatives in eq. (|2.8|) involve just the 
colour source pr = p{x~) in this outermost layer. 

To conclude, when integrating out the quantum modes in layers of r, one builds the 
colour source in layers of x~, with a one-to-one correspondence between the x~ coordinate 
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of a given layer in p and the rapidity r of the modes that have been integrated out to 
generate that layer. 

This correspondence is most simply formulated if one uses the space-time rapidity y, 
y = ln(x~/xo), Xo=l/P+, -cx)< y < oo , (2.9) 
to indicate the longitudinal coordinate of a field. We shall set, e.g., 
Py(x^) = x^p"(x^,x_l) for x~ = Xy = Xq e^, 
rfypy(x±) = J dx- p''ix-,x^), (2.10) 

and similarly for a, eq. ( |2.7] ), or any other field. Eq. ( p.3[ ) is then rewritten as : 

^.(.J^Pexp{../.y„;,..)r}. (2.U) 

Then, the previous discussion shows that the space-time rapidity y of a given layer 
in p is identical to the usual (momentum) rapidity of the fast gluons that have produced 
that layer. In particular, the colour source created by the quantum evolution up to r has 
support at space-time rapidities y < r (in agreement with the simple argument based on 
the uncertainty principle). This is important since, as we shall see, the r-dependence of 
the observables computed in the effective theory comes precisely from the upper limit on 
the longitudinal support of p. The Colour Glass evolves by expanding in y. 



3 The Gaussian approximation 

The exact solution to eq. ( ^.8| ) is certainly not a Gaussian. The kernel x[p] in this equation 
is non-linear in p to all orders, so the complete solution VF,-[p] must be non-linear as well, 
for any initial condition. (This non-linearity is manifest on the formal representation 
of the solution derived in Ref. [0.) Still, as mentioned in the Introduction, there are 
interesting situations in which a weight function which is a Gaussian in p may capture 
the relevant physics. Before describing such physical situations in more detail, in Sections 

4 and 5, let us first characterize the general structure of a Gaussian weight function, and 
show how to obtain such a Gaussian approximation for WV[p]; at least at a formal level, 
via mean field approximations to the non-linear RGE ( |2.8| ) (see also Ref. 0]). 

We start by displaying the most general Gaussian weight function which is consistent 
with gauge symmetry, and respects the correspondence between space-time rapidity and 
momentum rapidity explained in the previous section. In the covariant gauge, in which 
both the classical solution and the kernel x[p] are known explicitly, this Gaussian reads: 
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where 5^- [p] is a ^-functional enforcing that py = for any y > r, and Mr is a numerical 
factor, which ensures the right normalization: 



jDpWr[p] = l, (3.2) 

but is irrelevant for the computation of correlation functions. 

Note that the information about the longitudinal support of Py is included in the 
weight function. Thus, when computing observables by averaging over p with this weight 
function, like in eq. (|2.2|), the support of Py will be restricted to y < r, as it should. 

For physical space-time rapidities y < r, eq. (|3.1|) is indeed a Gaussian in p, with the 
following 2-point function : 

(Py(^±) Py{y±))r = S^'6{y - y')e{T - y) X,{x^, y^), (3.3) 

which is local in space-time rapidity, but generally non-local in the transverse coordinates. 
This is the most general non-local structure which is allowed by gauge symmetry, as we 
explain now: For Wr[p] to be gauge-invariant, the colour sources Py(a;_L) at different 
points (j,x±) must be connected by appropriate gauge links, or "Wilson lines". In the 
covariant gauge, in which eq. (|3.1| ) is explicitly written, the gauge potential has just a plus 
component, A'^ = aa, so the only non-trivial gauge links are those in the longitudinal 
direction, cf. eq. ( p.3|) . And, indeed, the kernel x[p] of the RGE is built with such Wilson 
lines 0, and thus is consistent with gauge symmetry, but also non-linear in p to all 



orders. For the gauge-invariance to be preserved in the Gaussian approximation, where 
no Wilson lines are permitted, WV[p] must be local in y, as shown in eq. ( p.l|) . 
The Gaussian weight function (|3.1|) obeys the following evolution equation: 



By comparison with eq. ( |2.(j| ), one sees that, at a formal level, the Gaussian approximation 
is tantamount to replacing the general kernel x[p\ iii the RGE by the p-independent, but 
r-dependent, quantity A,-. 

To determine the kernel A,- of the Gaussian, which is also the 2-point function in the 
present approximation, we shall require this quantity to obey the "mean field version" 
of the actual evolution equation satisfied by the complete 2-point function (p(l)p(2))^. 
Here, by "mean field version", we mean that the correlations which enter this equation 
are evaluated with the Gaussian weight function ( ^.1| ) itself. This gives a self-consistent 
approximation for the 2-point function. 

Since the general 2-point function (Py(x^) Py,{y±))r is not local in y, nor diagonal in 
colour, it is preferable to consider the equation satisfied by the following quantity: 

firix^, y^) ^ (p"(xx) p\yi))r^ (3.5) 
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where 



p\x^) ^ j dyp^ix^) (3.6) 

is the effective colour charge density in the transverse plane, as obtained after integrating 
over the longitudinal profile of the hadron. 

In the Gaussian approximation, eq. (|3.3| ) immediately implies: 

dr 



PrixL.yA.) = / c?yAy(x_L,l/±) , or A^(x_L,l/±) = ^(a;^,Z/±). (3.7) 



In general, the evolution equation for the two-point function (^]5|) can be obtained 
from its definition: 

{p\x^)p\y^))r ^ j Dp WM p\x^)p\y^), (3.8) 

by first taking a derivative with respect to r, then using eq. (p.8| ) for dWr/dr, and finally 
performing a couple of functional integrations by parts. This yields: 

|- {p\x^) p\y^)), = {x^\x^, y^))r + (a'^(x^) /(y^)). + (p'^(x^) a\y^))r , (3.9) 
where the quantity a°-{xj_) has been generated as the functional derivative of x[p\'- 

It is useful to recall here that, in the quantum calculation leading to eq. ( ^.3] ) 0], a[p\ 
represents the "virtual correction" (i.e., the sum of the self-energy and vertex corrections), 
and is complementary to the "real correction" represented by x[p]- The relation ( |3.1U| ) 



between these two quantities [gO], ^ is very important, since it ensures that the evolution 
equations for gauge-invariant observables derived from eq. ( p.8|) are infrared finite. 

In general, both x and a are non-linear in p, so eq. ( p.9| ) is not a closed equation — it 
couples the 2-point function to higher n-point functions of p — , but only the first in an 
infinite hierarchy of coupled evolution equations. A closed equation can be nevertheless 
obtained in the mean field approximation (MFA) in which the correlation functions in 
eq. ( |3.9| ) are evaluated with the Gaussian weight function (p.l|). Then, eq. (|3.9| ) reduces 



to the following equation for = dpr/dr (cf. eq. (|3.71)) 

Kix±,y±) = {xix±,y±) + a{x±)p{y±) + pix±) a{yA_))r , (3.11) 

where (x) = (x"'^)/(^c ^^'^■^ expectation value in the r.h.s. is computed with 

the weight function (|3.1|) , so it is a known functional of Ay. In general, this functional is 
very non-linear and also non-local, but it considerably simplifies in the limiting cases of 
interest (see Sect. 5 below). 



11 



To summarize, our MFA is defined by the weight function ( ^.1| ) with the kernel A,- 
determined self-consistently by eq. ( |3.11| ). This is more general than the MFA considered 
in Ref. [Q, which was an approximation on the evolution Hamiltonian (the functional 
differential operator in the r.h.s. of eq. 



H 



1 A r 1 A 

2 5pr ^^^^ 6pr 



(3.12) 



In eq. ( |3.12|) , the kernel x[p] the Hamiltonian is simply replaced by its average over p, 
so its functional derivative a[p] is implicitly neglected, and eq. (|3.11| ) reduces to 



(3.13) 



as is also obvious by comparing eqs. ( p. 4]) and (|3.12|) . The reason why eq. (|3.11|) appears 
to be more complete, is that, in its derivation, the equation for the 2-point function has 
been first generated with the full Hamiltonian, and the MFA has been implemented only 
afterwards. By contrast, eq. ( |3.13| ) has been generated directly with an approximate form 
of the Hamiltonian. 

As we shall discuss in Sect. 5, the difference between eqs. ( p.ll| ) and ( p.l3| ) has 
important consequences. In particular, eq. ( p.l3| ) is not good enough for our present 
purposes. Rather, our final approximation for the weight function, to be constructed in 
Sect. 6, will be based on eq. (p.ll|) . 



We conclude this section by displaying the expression for the dipole-hadron scattering 
amplitude (|2.2| ) in the Gaussian approximation. The Wilson lines in eq. ( |2.11|) involve the 
Coulomb field a'^{x±), which is linearly related to the colour charge density (cf. eq. (|2.4|) ), 
and therefore is itself a Gaussian random variable, with 2-point function (cf eq. (|3.3|)) : 



(3.14) 



where 



f 1 1 

ly{x±,y_L) = Jd'^z±d'^u±{x±\—^\z±)Xy{z±,u_L){u±\—^\y±). (3.15) 



From eqs. (|2.11| ) and (|3.14|) , a straightforward calculation yields [0, |TOl : 



ST{x±,y_L) = exp 



R 



dj 



ly{x_L,x±) + -fy{y_L,y±) - 2-fy{x±,y±) }. (3.16) 



The colour factor Cr is Cp = {N^ — l)/2Nc for a dipole made of a quark and an antiquark 
(this is the relevant case for 7*p deep inelastic scattering), and Ca = Nc for a gluonic 
dipole, which is interesting for the discussion of the RGE (since the Wilson lines which 
enter the kernel x[p] ^^re written in the adjoint representation, i.e., they refer to gluons). 

Eq. ( |3.16| ) simplifies if one assumes the hadron to be homogeneous in the transverse 
plane, within a disk of radius R, with R much larger than any transverse separation x±—y± 
of interest (so that one can neglect the edge effects). Then, the various 2-point functions 
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depend only upon the relative coordinate: e.g., \y{x±,y±) = \y{r±), with = x± — y±, 
and it is convenient to introduce the Fourier transform: 



Ay(A;j 



-ik I -r 



and similarly for the other functions. The relation ( pj.l5 ) simplifies to: 



while eq. ( p. 161) is finally rewritten as: 



1 - e 



ikj_ ■rj_ 



(3.17) 



(3.18) 



(3.19) 



4 The McLerran-Venugopalan model 

In this section, we interrupt our presentation of the quantum evolution and briefly review 
a simple Gaussian model for the weight function, originally proposed by McLerran and 
Venugopalan ^ to describe the gluon distribution of a large nucleus (A S> 1) (see also 
0, |, H, 1^, 0). There is no quantum evolution in this model (i.e., no dependence upon 
x), but this can be viewed as the initial condition, valid at some moderate value Xq, for 
the quantum evolution towards smaller values of x. 

In this model, the only colour sources are the A x Nc valence quarks, which are 
assumed to be uncorrelated, since they do not overlap with each other: Valence quarks 
which overlap in transverse projection belong typically to different nucleons, and thus 
have no overlap in the longitudinal direction, because of confinement. The distribution of 
these sources is therefore described by a Gaussian weight function, like in eq. ( p.l| ), but 
with a kernel which is local also in the transverse coordinates: 



Py{x±)p^{x^) \ 



(4.1) 



Physically, the kernel Xy{x±) represents the average colour charge squared of the valence 
quarks per unit space-time rapidity and per unit transverse area, at the given point (y, x±). 
It satisfies the sum rule: 



dy / d'^x± Xy{x±) 



(4.2) 



which follows since the total colour squared in the nucleus is {Q"'Q°') = g'^CpANc (recall 
that the colour charge squared of a single quark is gH°'t°- = g'^Cp), with 



dy / d'^x±p^{x±). 



(4.3) 
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The y-dependence of the kernel Xy{x±) remains unspecified in the MV model. This is not 
necessarily a problem since, when computed in this model, many interesting quantities 
(see, e.g., eq. ( |3.19D ) involve only the colour charge integrated over y (cf. eqs. ( |3.5| )-(|37?|)): 

/TO 
dyXy{x±). (4.4) 
-oo 

To get simple formulae, it is again convenient to assume transverse homogeneity within 
the nuclear disk of radius Ra = -Ro^^''^ (this assumption can be relaxed, and actually is. 



in numerical simulations |^). Then 



/^o(^) = ^ « A'"- (4-5) 



When the MV model is used to compute physical quantities, one usually finds loga- 
rithmic infrared divergences, which can be attributed to the fact that the Gaussian weight 
function ( [4.1|) is local in x±. Note that, when translated to momentum space, this locality 
means that the Fourier modes of Ay are independent of fc^ (cf. eq. ( |3.17| )) : 

\{xL,yA_) = 5^^\x^-y^)\ =^ \{kA.) = \ ■ (4.6) 

Consider then eq. (|3.19|) for the dipole scattering amplitude. In the MV model, it becomes: 

^o(rx) = exp f [l - ^^'"■^"] } ' (4-7) 

and the integral over in the exponent has a logarithmic infrared divergence indeedQ. 

Clearly, this divergence refiects the limitations of the assumptions underlying the 
MV model. In particular, the assumption that the colour sources are uncorrelated must 
fail for sufficiently large transverse separations. For instance, long-range correlations 
will certainly occur because of confinement: The three valence quarks within the same 
nucleon are bound in a colour singlet state, so the total colour charge ([4.3|), together with 
its higher multipolar moments, must vanish when measured over distances of the order 
of the nucleon size i?o, or larger. This in turn requires the correlator \y{k^_) to vanish at 
low momenta /c_l ^ I/-R0, in contradiction with eq. (|4.6| ). 

Based on this physical idea, it has been proposed to cure the infrared problem 
of the MV model by imposing colour neutrality on the nucleon level. This means that, 
effectively, the integral over k± in eq. ( |4.7| ) is cut off in the infrared at momenta of order 
I/-R0 ~ ^QCD- In fact, since the divergence is only logarithmic, then, to leading-log 
accuracy, the exact value of the cutoff is not important, neither are the details of the 



non-perturbative mechanism responsible for this screening. To evaluate eq. ( [4.7| ) to this 



"'The divergence is only logarithmic since, at small k±, 1 — e*'''-^ '"-!- « —ik± ■ r± + {l/2){k± ■ r±)^, and 
the term linear in k± vanishes after angular integration. This compensation of the leading IR divergence 
between 1 and e*'^^ ''-^ is ultimately a consequence of the fact that the external probe is a colourless dipole, 
with a gauge-invariant S'-matrix element. 
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accuracy, one can insert a sharp momentum cutoff equal to Aqcd, and also expand the 
exponential (since k± ■ rj_ <^ 1 for this dominant contribution). One obtains 0, [1^ : 

S„(rjc.exp|-ri^ln-^}, (48) 

with: 

QUA) ^ a,Cnfio{A) = . (4.9) 

Still, for the physical scales of interest (e.g., in relation with heavy ion collisions), the 
leading-log approximation above is not really satisfactory, and the complete, numerical 
results in Ref. turn out to be quite sensitive to the specific assumptions about the 
non-perturbative physics (i.e., upon the prescription used to implement colour neutrality). 

Thus, in this classical context, where the valence quarks are the only colour sources, 
one cannot avoid a rather strong sensitivity to the soft, non-perturbative physics. As 
we shall show in what follows, this sensitivity disappears once the effects of the quantum 
evolution are taken into account: At sufficiently large energies (or small enough values 
of x), the colour sources are predominantly gluons, and these gluons have long-range 
correlations associated with the non-linear (but still perturbative) physics of saturation, 
which ensure colour neutrality already over the short scale l/Qs ^ V^QCD- 



5 Colour neutrality from quantum evolution 

We now return to the quantum evolution described by eq. ( p.8|) , and examine the validity 
of the Gaussian approximations obtained in Sect. 3 via formal manipulations. To this 
aim, we need to consider the structure of the kernel x[p] some detail. In fact, for the 
present purposes, it is enough to know that x depends upon p via the Wilson lines ( |2.11| ), 
and that the strength of the field a in these Wilson lines is correlated with the gluon 
density in the system (see Sect. 7.1). This suggests different approximations at short 
and, respectively, large distances compared to the saturation length l/Qs- 

5.1 High k_i : The BFKL equation 

For short distances r± <^ 1/Qs{t), that is, for an external probe with a large transverse 
resolution = l/rj_ (e.g., a small dipole with transverse size r_|_), scattering is dominated 
by gluons with relatively large transverse momenta {k± ^ Qsi^)), and thus a low density. 
In this dilute regime, the colour field is weak, ga <^ 1, so one can expand the Wilson lines 
( p.ll| ) in powers of ga, and keep only the linear term in the expansion: 

V^(xj_) ^ 1 + ig dYay{xj_) = 1 + iga{x±). (5.1) 
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After this expansion, x[p] becomes quadratic in p(x_l) = / dj py{xi_) (cf. eq. 

X[p] ^ pICp, (5.2) 
which means that the corresponding RGE is still non-linear : 

That is, even in the dilute regime at high k±, the weight function W^t[p] is, strictly 
speaking, not a Gaussian. Still, as compared to the general RGE ( |2.8|) , the evolution 
generated by eq. ( ^.3|) exhibits an important simplification: it does not mix correlations 
(p(l)p(2) . . ■ p{n)) with different numbers n of colour sources 

In particular, eq. (|5.3|) provides a closed evolution equation for the 2-point function 
( ^.51 ) — i.e., eq. ( |3.9| ) with x[p] ^ pJCp and a[p] ^ K, p — , which is moreover linear 



When written in momentum space, this equation is recognized as the BFKL equation ||13||: 



- d p±^— r^/ir(P±) 



dr 71^ J — P± 



- ^ / , Pr{k^) . (5.4) 

The emergence of the BFKL equation in this regime, and in the present approximations, 
is as expected, since the quantity Prikx) plays also the role of the unintegrated gluon 
distribution (see Sect. 7.1). The two terms in the r.h.s. are the "real" and, respectively, 
"virtual" BFKL corrections, and originate from the terms involving x and, respectively. 



0", in the r.h.s. of eq. (|3.9| ). Since eq. (|5^ ) is local in rapidity, then, clearly, if Pr{ki_) 
satisfies this equation, so does also Xr{k±^) = dpr/dr, which is the kernel of the Gaussian. 

Let us now compare eq. ( [5.4| ) with the corresponding predictions of the mean field 
approximations introduced in Sect. 3, namely, eqs. ( |3.11| ) and (|3.13|) . 

In this low density regime, where x is quadratic in p, while a is linear, there is obviously 



no difference between the general equation (p.9|) and its mean field approximation ( p. 11 ). 



Thus, the BFKL equation (|^) is exactly reproduced by the linearized version of eq. (PTTTf) 



That is, although the general solution Wr[p] to eq. (|5.3|) is not a Gaussian, our mean 
field approximation is still able to encode the evolution of the 2-point function correctly, 
because this evolution is linear, and it does not mix the 2-point function with higher n- 
point functions (i.e., because the BFKL Hamiltonian in the r.h.s. of eq. (|5.3| ) is diagonal 
in the number of fields p). 

It is also important to notice that, within this Gaussian effective theory, the BFKL 
evolution of the kernel Xr of the Gaussian gets transmitted to other 2-point functions of 
physical interest, like the dipole-hadron scattering amplitude (p^). Specifically, in the 
linear approximation appropriate at high k^, eq. (|3.19| ) reduces to: 
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By taking a derivative with respect to r in this equation, and using eq. ( p.4| ) for A,- = 
d/j^r/dr, one can check, after some lengthy but straightforward manipulations, that the 
(hnearized) scattering amplitude obeys to: 

^■AC(rx) = as 7 ^-K^ (^r{z^) - iKir^)] , (5.6) 



which is the coordinate space form of the BFKL equation. This is further recognized 
as the linearized version of the Balitsky-Kovchegov (BK) equation, which is the general 
evolution equation for N'r{r±) that follows from the RGE ( |2.8|) 0, and also from other 



formalisms which focus on the evolution of the dipole wavefunction [^, |2^, |30|, q1 



Consider also the MFA based on eq. (|3.13 ) : Since this approximation misses the 
a term, it fails to reproduce the "virtual" BFKL term, and thus the complete BFKL 
equation ( |5.4| ). This is a serious failure: as well known, the "virtual" term cancels the 
infrared divergence of the "real" term at p± = k±, and thus plays an important role in 
the BFKL dynamics. 

Still, there exists a kinematical limit in which the "virtual" term becomes unimpor- 
tant: this is high-/c^ regime where one retains only terms enhanced by both energy logs 
ln(l/a;) and transverse momentum logs ln(fc^/A^) ("double logarithmic approximation", 
or DLA). This limit is common to the BFKL and DGLAP |^ equations, and is 
formally obtained by assuming p±_ -C k±_ in eq. ( |5.4| ). In this limit, one expects the MFA 
based on eq. ( 3.13| ) to become appropriate. This is confirmed by the analysis in Ref. @] 
which shows that, at very high k±, the self-consistency equation ( p.l3| ) reduces to: 

K{k^) = ^ r dj f'^Xyip^), (5.7) 



J-oc J P± 

or, equivalently (cf. eq. (|3.71 )), 

dfirik^) asNc [^^^ dp 



, 2 ^^Ap±), (5.8) 

which is recognized indeed as the DLA evolution equation (compare to eq. ( |5.4| )). 

It is important to keep in mind, however, that DLA is not the right approximation to 
study the approach towards saturation. Indeed, there is a large kinematical gap between 



the validity range for this approximation and the saturation line k± = Qs(t) [T^, and, 
within this gap, one can rely only on the full BFKL equation. This is why in what follows 
we shall focus exclusively (in this high momentum regime) on the BFKL dynamics. 



5.2 Low k± : Saturation and colour neutrality 

At low transverse momenta kj_ < Qg, the gluons are saturated, that is, they have large 
occupation numbers dN / dT(Pk± ~ l/a<i, and a radial momentum distribution dN/dTdk± 
which is peaked at the relatively hard scale /cj_ ~ Qs{t) (see Sect. 7.1). Thus, in the 
corresponding CGC description, the colour fields have large amphtude^ a ~ 1/g, and 

^More precisely, we mean here the fields integrated over all rapidities, i.e., a{x±) ~ Jdyay{x±). 
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carry typical momenta k± ~ Qg. Accordingly, the Wilson lines ( |2.11 ) — which are 



complex exponentials built with these fields — oscillate around zero over a characteristic 
distance ~ 1/Qs{t) in the transverse plane. This implies that Wilson lines which are 
separated by large distances ^ 1/Qs{t) are necessarily uncorrected (since their relative 
phases are random). Thus, when studying the dynamics over large transverse separations 
^ 1/ Qsi^), it should be a good approximation to neglect the correlations of the Wilson 
lines (or, more generally, to treat them as small quantities). This is the "random phase 
approximation" (RPA) introduced in Ref. [§, 0. 

In this approximation, the RGE ( p.8|) simplifies drastically : After neglecting the 
Wilson lines, the kernel x becomes independent of p, and the corresponding solution Wr[p] 
is truly a Gaussian (to the accuracy of the RPA). This is quite remarkable, given that, 
at saturation, we are in a strong field regime, where the dynamics is fully non-linear. 
Specifically, one finds p 

X^\x^,y^) ^ 5'^'-{x^\-Vl\y^) (RPA) (5.9) 

TT 

which gives the following RGE : 

dWr[p] _ 1 fd^k^ kl 5^Wr[p] 



dr 2 J (27r)2 tt 6p'^{k^)6p'^{ 



(5.10) 



where we recognize a Coulomb-like Hamiltonian in the r.h.s. A brief comparison with 
eq. (|3.4|) allows us to identify the charge-charge correlator at low k± : 



\y{k^) = -kl, (A;^«g.(r)). (5.11) 

71 



Clearly, both mean field approximations discussed in Sect. 3 reproduce eqs. (|5.9| )- (|5.11 
in the low-k± limit. 

Note that this low-fcj, distribution is homogeneous in all the (longitudinal and trans- 
verse) coordinates. That is, in the transverse plane, it is only a function of the relative 
coordinate x_i — y±, and in the longitudinal direction, it is independent of the space-time 
rapidity y. One should however keep in mind that, for given kj_ -C QsiT), eq. ( |5.11| ) 



applies only for y in the interval Ts{k±) < y < r, with Ts{k±) being the rapidity at which 
the saturation scale becomes equal to the momentum k± of interest: 

g2(r) = kl for r = r,(A;^). (5.12) 

It follows that the integrated quantity (cf. eq. (p^)) : 

ftr{k±) = [ dy^ =(r- Ts{k^)) ^ , (k^ « a(r)), (5.13) 

sat y 77 V / TT 

which measures the density of saturated colour sources (with given k±) in the transverse 
plane, grows only linearly with r, that is, logarithmically with the energy. This should 
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be contrasted with the rapid increase of the corresponding quantity at k±_ ^ Qsij)i 
which evolves according to the BFKL equation (|5.4]) , and thus rises exponentially with r. 
Eq. ( |5.13| ) shows that, at low momenta k±_ <C Qs(t), the colour charge density saturates, 



because of the strong non-linear effects in the quantum evolution [R|. As we shall discuss 



in Sect. 7.1, this further implies the saturation of the gluon distribution. In Ref. [^, the 



result in eq. (5.13) has been reobtained from a study of the BK equation. 



According to eq. (|5.11| ), the only remaining correlations at saturation are those in 



the transverse separation r_L between the colour sources. Importantly, these are such as 
to ensure colour neutrality ||, Indeed, the fact that Ay(fc^) oc k\ when k\ ^ Q 
implies {Q^Q"") = 0, where Q"" is the total colour charge ( |4.3| ) of the saturated gluons: 



(Q'^Q'^). = J dyXyik^ = 0) = 0. (5.14) 

Since the distribution of p at small k± is a Gaussian, and both the 1-point function 
{Q"'), and the 2-point function of vanish when computed with this Gaussian, 

it follows that = 0, which physically means that the saturated gluons are globally 
colour neutral. In fact, since the Wilson lines average to zero over a finite distance 
Ax± ^ 1/Qs{t), it follows that colour neutrality is achieved already over a transverse scale 
of the order of the saturation length (this can be explicitly verified by using eq. ( p.ll| )): 

[dy [ c/Vpy(xx) = 0, (5.15) 
J Jasj_ 

where ASj_ is, e.g., a disk of radius R > 1/Qs{t). 

When k± ^ 0, the spectrum ( [5. 11] ) vanishes sufficiently fast to ensure the infrared 
finiteness of the dipole scattering amplitude ( p.l9|) , and, more generally, of all the quan- 
tities which were, at most, linearly infrared divergent in the MV model. Thus, we expect 
that all the observables computed within the CGC effective theory come out infrared 
finite, and thus are insensitive to the non-perturbative physics at momenta k± ^ Aqcd- 
In fact, because of saturation, the spectrum starts to soften already at momenta k± > Qg, 
where the linear, BFKL, evolution applies. This happens in the window for "geometric 
scaling" [13], to be discussed in Sect. 6. 



5.3 Brief summary of the quantum evolution 

To summarize, if we compare to the initial condition of the MV type — namely, the 
local Gaussian ([4.1|), which describes a system of independent colour sources — , the 
general effect of the quantum evolution is to introduce correlations in the transverse and 
longitudinal directions (i.e., non-localities in r± and y), and non-linearities in p (i.e., 
higher ra-point correlations among the colour sources). However, the role and the relative 
importance of the various non-local effects, and of the higher ra-point functions, may 
differ from one kinematical regime to another. 
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At transverse momenta far above Qs-, the longitudinal correlations are irrelevant (they 
are always integrated out, as in eqs. ( |3.6|) and (pTT|)), and the various n-point correlations 
evolve independently from each other. In particular, the 2-point function evolves accord- 
ing to the BFKL equation (p.4|), and so does the dipole-hadron scattering amplitude. 

In the opposite limit at /c± ^ Qs, the RGE becomes linear (cf. eq. (|5.1CI| )), and the 
weight function is strictly Gaussian, and therefore local in y. The only correlations which 
persist in this regime are the transverse correlations reflecting the colour neutrality of the 
saturated gluons. These are sufficient to ensure that gauge-invariant observables come 
out infrared finite when computed in the effective theory. 

In the intermediate regime at fcj. ~ Qs, the evolution is fully non-linear — correlation 
functions of all orders get mixed in the evolution equations — , and therefore very non- 
local in y. Clearly, there is no hope to describe the complex physics in this transition 
regime via a Gaussian approximation. However, we expect this physics not to be essential 
for physical quantities for which the external resolution scale is either well above, or well 
below, the saturation scale. As we shall verify on the example of the dipole-hadron 
scattering amplitude in Sect. 7.2, such quantities receive most of their contributions from 
either the hard modes with Qs, or the saturated gluons with k±_ <^Qs- 

To easily compute such observables, it would be advantageous to dispose of a Gaussian 
weight function which interpolates smoothly between the known behaviours at high and, 
respectively, low transverse momenta. For instance, the complete solution Ar(xj_, to 
the evolution equation in the MFA, eq. ( p. Ill ), would provide such an interpolation, since 
this is necessarily a smooth function and has the right limiting behaviours. But eq. ( |3.11|) 
is a complicated, non-linear, integral equation, whose general solution can be obtained, 
at most, numerically. For practical calculations, it would be preferable to have a simple 
analytic interpolation, which satisfies the BFKL equation (|5.4| ) at high /c^, and reduces to 
eq. (|5.11| ) at low A;^, with the transition between the two regimes taking place at k± ~ Qs- 
Such an interpolation will be constructed in the next section. 



6 The interpolation 

So far, we have found the following limiting behaviours for the charge-charge correlator 
( p.3|) at high and, respectively, low transverse momenta: 

^ ' \ ^kl, for ks_ < Qsir) ^ ' 

Here, and from now on, we reserve the notation At-(/c±) for the solution to the BFKL 
equation, and use Xrik^^) to denote the kernel of the Gaussian weight function (|3.1|) . 

Clearly, the following function provides a smooth interpolation between the limiting 
behaviours shown above: 
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There is, however, one more constraint that eq. ( |6.2| ) must satisfy in order to be a physi- 
cally acceptable interpolation: the change from the high-momentum regime to the low- 
momentum one must happen at the saturation scale. For the function in eq. ( |6.2| ), this 
condition is satisfied provided the two terms in the denominator become of similar mag- 
nitude when k± = Qsi^) '■ 

-Ql{r) ^ \r{k^ = Qs{r)). (6.3) 

TT 

Note that this is precisely the condition that the limiting behaviours shown in eq. (|6.1|) 
match quasi-continuously with each other at the saturation scale. Clearly, whether this 
condition is fulfilled or not, depends upon our definition of the saturation scale. Usually, 
Qs{t) is defined through a "saturation criterion", imposed on the gluon "packing factor" 



in the transverse plane [14, 15, 16, 17 1, or, equivalently, on the dipole-hadron scattering 



amplitude |^ |2^, |, |TB|, |T^, ^, For instance, in terms of A/'r(rx), the saturation 
criterion reads (k is a number of order unity) 

Mr{r^) = K for ri = l/Q2(r), (6.4) 

which is the same as the condition that the dipole of transverse size equal to the saturation 
length is completely absorbed ("blackness"). 

We shall argue below that, in fact, eq. ( |6.3| ) is just another definition of the saturation 
scale, equivalent to that in eq. (|6^ ). In other terms, the condition (|6.3| ) is satisfied at the 
right saturation scale hy definition, so eq. ( |6.2| ) is a physically acceptable interpolation 
indeed. Of course, such an interpolation is not unique, but at the present level of accuracy 
we cannot distinguish between different interpolations which are all physically acceptable 
according to the previous criteria. In what follows, we shall adopt the interpolation (|6.2|) 
because of its simplicity. 

At a mathematical level, the equivalence between the two formulations of the sat- 
uration condition — in terms of the scattering amplitude, eq. (|6.4|) , or in terms of the 
charge-charge correlator, eq. ( |6.3D — is almost obvious: Both A/'r('"±) and At-(A;_l) satisfy 
the BFKL equation and obey the same saturation criterion, namely, the relevant "scaling 
function" (A/'r(r_|_) in one case, and \T{ki_) /k\ in the other) is constant on the saturation 
line = Qs{t). Thus, it is a priori clear that eq. (|6.3| ) will provide the correct estimate 
for Qs{t). Still, it is instructive to briefly develop the consequences of this equation here, 
in order not only to recover the known result for the saturation scale [|17|, |1^, but also to 
establish a simple scaling approximation for At-(A;_l), which will greatly simplify the use 
of eq. ( |6.2| ) in practice. 

For convenience, in what follows we shall mostly consider the solution to the BFKL 
equation ( |^ ) for the case of a fixed coupling a^, and only briefly indicate the changes 
which occur when a running coupling as{Ql) is used instead. Also, we shall need only the 
approximate form of the solution, which is obtained via a saddle point approximation in 



the Mellin transform of At-(/c_l) (see, e.g., |jT8[ for details): 

2- 



Xr{k^) ^ /(r)y/gg^exp|^«,r-^^(^ln^^ |. (6.5) 
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In this equation, /(r) is a slowly varying function that will be ignored in what follows 
(since its contribution to In^X^/k"^) is comparatively small), Qq is a scale of order Aqcd 
introduced by the initial conditions, ds = C(sNc/n, and u and /3 are pure numbers: u = 
41n2 ^ 2.77, /3 = 28C(3) ^ 33.67. 

Eq. ( |6.5D is strictly correct for k± ^ Qs{t), where the colour charge density is low 
[Xr/ k"^ <^ 1), and BFKL is a right approximation. But this expression can be extrapolated 
down to k±_ ~ Qs{t) (since BFKL remains marginally valid down to the saturation line), 
and then inserted in the saturation condition ( |6.3D , to obtain an equation for Qs{t) : 




exp \ LuasT - _ ( In — I ) \ ^ I. (6.6) 



This has the following solution : 

Ql{r) = Q^e'="^^ c = [-/3+y^pTM]/2 = 4.84... (6.7) 
which is the correct result IIH], as anticipated. Note that the factor in front of the 



exponential in eq. ( |6.7|) is not under control in the present approximations. A more refined 
treatment shows that this prefactor has actually a weak dependence upon r. 

It is further interesting to study the behaviour of Ar(fc±) for momenta close to the 
saturation scale, where a simple approximation to eq. ( |6.5|) can be obtained. To this aim, 
let us replace Ql by Q1{t) as the reference scale in eq. (|6.5|), by writing: 

A simple calculation yields then: 

Kik.) ^ i*i exp { J-^ - ^ (.„ J-^)'^ }. (O) 

where 7 = 1/2 + c/ 13 0.64, and the (a priori, uncontrolled) prefactor has been chosen 
in such a way that the matching condition ( |6.3| ) is fulfilled exactly^. 

The last equation suggests a remarkable simplification: Assume that k± is sufficiently 
close to Qsij) (although still above it) for hiikPi/Qlij)) <^ a^r. Then, the second term 
in the exponent can be neglected compared to the first one, and eq. (|6.9|) reduces to 



K{k^) ^ - kl(^] , (6.10) 



which shows geometric scaling : the dimensionless function A^(A;^)/fc^ depends upon 
its two variables kj_ and r only via the combination k^/Q^i^r). Eq. ( |6.10| ) is a correct 
approximation for momenta fc^ within the following "scaling window" [III 



Qlir) « kl «^. (6.11) 



^Tliis is, of course, just a matter of convenience; an approximate matching, as shown in eq. (|6.3|) , 
would be enough for the present purposes. 
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The scaling form ( |6.10| ) is similar to the one derived in Ref. [|T^ for the dipole-hadron 
scattering amplitude, which reads (for within the scaling window (|6.11|) ): 



(6.12) 



with the same exponent 7 as in eq. ( |6.10|) . This exponent (or, rather, the difference 1 — 7) 
is sometimes referred to as an "anomalous dimension" |T^ , In fact, as we shall 

check in Sect. 7.2, the similitude between eqs. ( |6.1CI| ) and ( |6.12D is consistent with the 
linear relation (|5.5|) between these two quantities, valid for small r^. 

At this point, we should recall that geometric scaling has been first identified in 
the experimental data, as a regularity of the HERA data for 7*-proton deep inelastic 
scattering [^. More recently, this notion has found phenomenological applications also 
to electron-nucleus DIS [^], and to heavy ion coUisions . 
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Figure 1: Momentum dependence of Ay (A;_l) = {'n'/kj_)Xy{k±), plotted as a function of 
k±/Qs{T) for two values of 7. Solid (red) line: 7 = 0.64, and dashed (green) line 7 = 1. 



The previous considerations show that the following scaling Ansatz for the interpo- 
lating kernel (|6.2|) should be a good approximation for all momenta k± ^ Qlij) / Kqcd'- 



Ay (A; J 



^ 1 + 



^(y)V 



^(y) 



(6.13) 



We have reintroduced here the 6'-function to emphasize that the charge-charge correlator 
vanishes for space-time rapidities y > r. The quality of this interpolation will be tested 
in the next section, where we shall use eq. (|6.13|) to compute the dipole-hadron scattering 
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Figure 2: Longitudinal structure of \^y'^{k_[_) = {'n'/kj_)\y{k±). This is plotted as a 
function of y for two transverse momenta: k^^^ = 3Q'^{t) (green line, dashed), and 
k^l"" = 0.02QI{t) (red line, sohd). We have used 7 = 0.64, c = 4.84, a, = 0.2, and 
T = 10. The corresponding separation rapidities are (fcl^'s^) ~ 11 and Tsik^"^) ~ 5. 




amplitude and the unintegrated gluon distribution, and then we shall compare the results 
to known general properties. 

In such applications, it will be convenient to treat the "anomalous dimension" 7, 
and also the exponent c in eq. ( |6.7| ) for the saturation scale, as free parameters. For 
instance, eq. ( |6.13|) with 7 = 1 provides a simple generalization of the MV model (in the 
sense that the transverse correlations die out at high k±), but which includes the correct 
infrared behaviour due to saturation (thus avoiding the infrared problem of the original 
MV model), and also some of the quantum evolution, via the rapidity dependence of the 
saturation scale. 

Moreover, eq. ( |6.13|) with 7 = 1 can be used as a rough approximation for the general 
kernel ( p.2[ ) in the DLA regimeQ at very high transverse momenta k± ^ Q1{t) / Kqcd- In- 
deed, in that regime, the solution to the BFKL equation can be replaced by its DLA 
approximation, which has only a weak dependence upon fc^. This suggests that, in 
applications to the phenomenology, one can always use the simplified kernel ( |6.13|) , 
but let the anomalous dimension 7 be weakly dependent upon fc_L, in such a way that 
7 K, 7BFKL — 0.64 within the scaling window ( |6.11|) , and 7 ?a 1 for much larger momenta. 
More generally, one can imagine extracting both 7 and c from fits to the experimental 

^Note that this is the regime where the dipole-hadron scattering shows "colour transparency" i.e., 
the corresponding scattering amphtude behaves Hke Mr{r±_) ~ a^r^ xG{x, l/r'\_)/'KR^ , where xG{x, Q'^) 
is the gluon distribution to be discussed in Sect. 7.1. 
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data. In the subsequent applications, we shall give results for both 7 = 1 and 7 = 0.64. 

The quantity Xy{k±) is represented in Fig. |I] as a function of kj_ for fixed y, and in 
Fig. as a function of y for two values of k± (one above Qs{t), the other one below it). 
For convenience, we have plotted the rescaled function X^y'^{k±) = {'7i/k'^)Xy{k±). Note 
that the y dependence of Xy{k±) represents the longitudinal profile of the colour charge 
distribution in the hadron, for modes of given k_i_. (Recall that Xy{k±) is the average colour 
charge squared per unit space-time rapidity per unit transverse area, for a given k±.) If 
k± > Qs{t), the colour charge density increases exponentially with y all the way up to the 
edge of the hadron (located at y = r). li k± < Qs{t), there is an exponential increase up 
to the intermediate rapidity y = Ts{k±), where the gluon modes with momentum k± start 
to be saturated (cf. eq. ( ^.12| )); then, for rapidities Ts{k_i_) < j < r, the charge density 
remains nearly constant. This behaviour is manifest on Fig. |^. 

To conclude this section, let us indicate how the previous results change when, instead 
of a fixed coupling a.^, one rather uses a coupling which is running with the saturation 
momentum: — ^ ^^(Q^), with a<j(Q^) = bo / In^Q"^ / Aqq^^) . As shown in Refs. [1^, p^ , 
the only change refers to the functional form of Q1{t), which now becomes Q1{t) = 

^QCD^^^'''''^*''^''"'^"'' 5 where c is the same number as in eq. ( p.7|) , and tq is a constant. All 
the other results related to scaling — the value of the anomalous dimension 7 and the 
momentum range ( |6.11| ) in which the scaling holds — remain unchanged, except for the 
expression of the saturation scale which enters these results. In the applications to be 
considered in the next section, we shall restrict ourselves to the case of a fixed coupling. 



7 Some applications 

In this section, we shall consider a few applications of the Gaussian effective theory 
with kernel ( [6.131 ). First, we shall discuss the unintegrated gluon distribution, for which 
we shall derive a simple analytic formula which interpolates between saturation at low 
momenta and BFKL evolution at high momenta, and which can be easily used in further 
applications. In this context, we shall briefly recall how saturation arises in the effective 



theory for the CGC |]I0[, and emphasize the difference in this respect between the classical 
MV model ^ and the full theory including quantum evolution Then, we shall 
use eq. (|3.19|) to compute the dipole-hadron scattering amplitude and study its various 



limits. The resulting expression has the right scaling behaviour at short distances, and the 
correct approach towards "blackness" at large distances, and can be easily implemented 
in phenomenological studies of deep inelastic scattering. 

7.1 The unintegrated gluon distribution 



From eqs. (|6.13| ) and (|6.7|) , it is easy to derive the following approximation for the quantity 



I^T{k±) which characterizes the distribution of the two-dimensional colour charge density 
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p°(ix) = Idyp°(x^) (cf. eqs. (g^-Q) : 

This quantity is interesting for, at least, two reasons: 

First, it determines the probabihty distribution for the colour charge density in the 
transverse plane (i.e., the three-dimensional distribution integrated over y) : 

Z,|/(.J|=A/-.exp|-i/ P'iM^) \ 

This is the relevant weight function for applications which do not consider explicitly the 
longitudinal structure of the colour source (like the numerical simulations in Refs. ||3^ ). 



Second, eq. ( [7.1D provides a good approximation for the unintegrated gluon distribu- 
tion in a wide range of transverse momenta. More precisely, /iT(^±)/^i coincides with the 
gluon distribution computed in the CGC effective theory at both high and low momenta 
compared to Qs(r), and it provides a smooth interpolation between these regimes. 

To show this, we start with the definition of the gluon distribution as the Fock-space 
gluon density in an appropriate gauge. Consider the canonical quantization of the Yang- 
Mills field theory in the light cone gauge = ||4^. The gluon distribution xG{x,Q'^) 
(= the number of gluons of transverse size Aa;^ ~ l/Q per unit rapidity) is then obtained 
as P, [ig] (with k+ = xP+, r = ln(l/a;)) : 



xG{x,Q') = jd'k^Q{Q'-kl)k+^^^. (7.3) 

The quantity d^N/d^k, which in eq. ( |7.3| ) plays the role of the unintegrated gluon distri- 
bution, is the Fock space gluon density, which is defined in terms of the gluon creation and 
annihilation operators in the standard way. This can be related to the following 2-point 
function of the LC-gauge "electric field" : 

ci^A^ 1 • 

(F^+(x+, fc^)Fr (x+, -k^)),, (7.4) 



which is in turn recognized as the LC-gauge expression of a gauge-invariant correlator, 
which in other gauges would involve also Wilson lines in the longitudinal direction ^ [ly] . 

In fact, in the effective theory for the CGC, one exploits this gauge freedom to rewrite 
the correlation function in eq. (7^) in terms of the gauge fields in the covariant gauge. 



for which we know the weight function. Namely, we use eq. ( p.5|) to write: 

{K^{S)F:^{y))r = {iUld^a%iUld^a%)r (7.5) 

with X = {x^,x±), and W{x) and aa{x) given by eqs. ( PTB| ) and (|2.7|), respectively. The 
Wilson lines in the r.h.s. of eq. ( [7.5|) encode all the non-linearities coming from the 
classical field equations (|2.1| ). 
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Within the Gaussian effective theory with weight function ( ^.1| ), this correlation func- 
tion can be exphcitly computed, with the following result M, ffH] : 



r dj I d'r^e''--^- S,{r^) (-Vi7y(r^)) , (7.6) 



where 'jy{r±) is the 2-point function of the field a, eq. (|3.15|) , and Sy{r±) is the 2-point 
function of the Wilson lines in the adjoint representation, i.e., eq. ( p.l9| ) with Cr = Nc- 
To simplify notations, it is useful to define: 

^ - 47r3 1 d^N 

a quantity that we shall refer to as the "unintegrated gluon distribution" in what follows. 
Up to the factor Att^, this is the number of gluons of each colour per unit rapidity per 
unit of transverse phase-space. 

Clearly, even with the interpolation ( |6.13D for the kernel of the Gaussian, which implies 



(cf. eq. ( |3.18| )) it is still not easy to perform the remaining integrations in eq. ( |7.6| ) 
analytically. But simpler formulae can be obtained in the limiting regimes where k± is 
either well above, or well below, the saturation scale. 

Consider the high-Zc^ regime first. Then, the integral in eq. ( |7.6| ) is dominated by 
small r± <C 1/Qs{t), so that 5'y(rj_) ^ 1, and the r.h.s. of eq. (|7^) reduces to 

r dj [ dV^e^^-- (-Vi7y(rx)) = r dj khyik^) = (7.9) 

J —oo J J —oo "'_L 

where we have also used eqs. (|3.18|) and ( p.7|) . This gives: 

V;.(fc^) ~ for A;x»gs(r), (7.10) 

with fir{kj_) obtained as the solution to the BFKL equation ( |5.4| ). Eq. ( |7. 10| ) is the 
standard BFKL approximation for the unintegrated gluon distribution. 

The \ow-k± case is more subtle, and we shall construct our argument in several steps: 
(I) We shall see eventually that, when k± ^ Qs{t), the dominant contribution to 
eq. ( [7.61) comes from space-time rapidities in the range Ts{k±) < j < t, with Ts{k±) 
defined in eq. ( ^.12] ). In order to check this, we first evaluate the contribution of the 
lower rapidities, y < Ts{k±). In this range, Qsiy) <^ k±, and the exponential factor e'"''-^'^-^ 
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selects distances r± ^ l/k± -C 1/Qs{j), for which Sy{r±) ~ 1 to a very good accuracy. 
Thus, this fow-y contribution is proportional to (cf. eq. ( [7.9|) ) 



1 



T=Ts{k^) as 



where the last estimate comes from using eq. ( |7. 1| ) . The contribution ( |7.11| ) is a priori 
large, since of order l/a^, but nevertheless subdominant, as we shall see. 
(II) We now consider the remaining contribution, namely: 

r dy I dh^e'^--^- S,{r^) (-Vi7y(r±)) , (7.12) 

where, for a typical y, k± ^ Qsij) Qs{t). Since e^'^-^'^^ selects relatively large distances 
r± ^ l/Qsij), we must consider the behaviour of the integrand for such distances. In 
momentum space, we have (cf. eqs. ( p.l8| ) and ( |5.11| ), or directly from eq. (|7.8|)) : 

khyik^) ^ l/vr, for ki_ < g,(y), (7.13) 

showing that, in coordinate space, — V5_7y(r^) = irS^'^^r^). Since eq. ( [7.13| ) is just a low- 
momentum approximation, this simply means that the function — V^7y(rx) is localized 
at small distances r± ^ 1/Qs(y), with integrated weight 

J d\^ (-Vi7y(r±)) = l/vr. (7.14) 

Thus, in reality, the integration in eq. (|7.12|) is restricted to values ^ 1/Qs(y), where we 
cannot compute the integrand accurately (since our approximations break down around 
Qs), but where Sy{r^_) 5 < 1 is essentially a number of order one. That is, we shall 
approximate: 

r dy I rfV^e^^--- Sy{r^) (-Vl7yK)) ^ - T dy = - (r - r,(A;^)) ,(7.15) 

where the number 5 < 1 cannot be computed in the present approximations. Since, 
moreover (cf. eqs. ( |6.7D and ( p.l2| )): 



where, by assumption, ln((5^(r)//c^) ^ 1, it is clear that the contribution ( [7.15|) is larger 
than that in eq. (|7.11| ), and thus is the dominant contribution, as anticipated. 
To conclude, at low momenta k_[_ <^ <5s(t), we have found: 

i^,{k^) ~ 5 ^ In , for k^ < g,(r) , (7.17) 
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which, up to the factor S, is the same as fiT{k±)/k'j_, with fir{kj_) given by its low- 
momentum approximation (|5.13|) . 



Eq. ( [7.17| ) exhibits gluon saturation. To see this, compare eqs. (|7.1CI| ) and (|7.17| ): At 



high momenta, (pT{k±) shows an exponential increase with r (according to the BFKL 
equation), and also a rapid increase with decreasing k± (as a power of l/k±), which, 
if extrapolated at low momenta, would generate infrared problems. But for k± below 
Qs{t), (pT{k±) is rather described by eq. (|7.17|) , which shows only a /^near increase with 
r, that is, logarithmic with 1/x — this is small-x saturation — , and also a logarithmic 
increase with l/k±_ — this is k±- saturation. The fact that both types of saturation 
occur simultaneously is consistent with the notion of "geometric scaling" : At saturation, 
there is only one intrinsic scale in the problem, the saturation scale, so ipr{k±) can be 
only a function of the ratio Q^t) / k^. Thus, a logarithmic dependence upon k^ entails 
necessarily a similar behaviour in r oc In QI{t). As we have seen in Sect. 6, the scaling 
property characteristic of saturation is preserved by the BFKL evolution until relatively 
large values of k±. 

At this point, let us also emphasize a crucial difference between the saturation mecha- 
nisms in the MV model and in the full effective theory, with quantum evolution included. 
In the MV model, /cj_-saturation occurs as well 0, but only as a consequence of the 
non- linear effects in the classical field equations ( p.l| ). That is, this classical saturation 



relies crucially on the presence of the Wilson lines in eq. (|7.5|) , or of their 2-point function 
S'y(r^) in eq. (^). 

On the other hand, these Wilson lines played almost no role in the quantum calculation 
leading to eq. ( [7. 171 ) (they are just responsible for the overall factor S). In this case, 
the saturation is rather the consequence of the long-range correlations among the colour 



sources, as encoded in the \ow-k± behaviour of the function Urik^) (cf. eqs. (|5.13|) and 
( [7.16| )). These correlations ensure the global neutrality of the colour sources at the scale 
1/Qs{t), and this in turn implies the saturation of the radiated gluons independently of 
the non-linear effects in the classical field dynamics. The relevant non-linear effects are 
rather those in the quantum evolution, i.e., in the RGE ( p. 81) , which generate the scale 
Qs{t) in the first place, as the scale for colour neutrality. 

These considerations explain why there is no contradiction, in this quantum context, 
between having an expression for the unintegrated gluon distribution which is linear in 
l^riki) both at high momenta and at low momenta, as we have found before, and which 
nevertheless shows saturation. This further implies that it is enough to dispose of a good 
approximation for the 2-point function /ir(/^±), namely eq. ( |7.1| ), to immediately deduce 
an equally good approximation for the unintegrated gluon distribution, that is. 



12/ 



T 



^^^In^lH-^^^J J. (7.18) 

At low momenta, k± -C <5s(t), this clearly reproduces the behaviour in eq. (|7.17|) (except 
for the factor 6, which is unknown anyway, and will be ignored in what follows). At high 
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Figure 3: <^r(fc±) as a function of k±/Qs{T) in linear scale (red line, solid), compared 
with its high momentum approximation (QHt) / k\)'' (green line, dashed). We have used 
as = 0.1 and 7 = 0.64. 



momenta, fc_L ^ Qs(t), 

1 ^Qlir) 



n^ycas 



(7.19) 



which is the geometric scaling version of the BFKL expression in eq. (|7.10|) , and therefore 
holds, strictly speaking, only within the scaling window (|6.11| ). 

According to ( [7.18| ), the saturation condition can be also formulated in terms of ipr '■ 
Qs{t) is the scale at which the unintegrated gluon distribution becomes of order l/a^: 

<^.(fc^) ~ — for k^^Qsir). (7.20) 



a.9 



The A;_L-dependence of the function ( |7.18| ) is illustrated in Figs. |^, | and |^, by using 
both linear and logarithmic scales, so that the various regimes alluded to before become 
manifest. In particular, in Fig. |^, we have indicated also the r-dependence of the result, 
by plotting iprikx) as a function of k±_ for several values of r. 



7.2 The dipole— hadron scattering amplitude 



After using eq. (|7.1| ), the ^-matrix element for dipole-hadron scattering, eq. ( p.l9|) , be- 
comes [for a gluonic dipole, i.e., with Cr = N^] : 



exp < - 



47r [(fk\ 1 



7c j (27r)^ 



In 



1 + 



(7.21) 
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Interpolation with y = 0.64 
Interpolation with y = 1 .00 
High mom. approx. 7= 0.64 
High mom. approx. y = 1 .00 

Low mom. approx. 




Figure 4: (^r(fc±) as a function of k^/Qsir) in log-log scale for two values of 7 : 7 = 
0.64 (red line, thick solid), and 7 = 1 (green line, thick dashed). Comparison is made 
with l/k^^^ 7 = 0-64 (magenta line, thin dotted), l/k^ (blue line, thin dashed), and 
(l/cis) ln(Q^(r)/A;^) (cyan line, thin dot-dashed). 



Clearly, for any 7 > 0, the integral in the exponent is both ultraviolet and infrared 
finite, so, by dimensional analysis, its result must be a scaling function (i.e., a function of 
^i.Qs{'t))- In what follows, we shall study the asymptotic behaviour of this function for 
both r^Qsir) <^ 1 (small dipole) and r^Qgij) ^ 1 (large dipole). This is interesting for 
several reasons: i) It allows us to identify the dominant scattering mechanisms in these 
limits, a ) It illustrates the role of saturation in providing a smoother infrared behaviour, 
and thus infrared finite results (this will be especially clear in the comparison with the 
MV model for 7 = 1). in) For a small dipole and 7 < 1, we shall verify explicitly that 
the anomalous dimension 7 is correctly transmitted from the kernel of the Gaussian to 
the scattering amplitude (cf. eq. ( |6.12| )). iv) For a large dipole, we shall recognize the 



behaviour expected from previous studies of the BK equation [^, |TT|, and of the RGE 
( p.8[ ) H, |1^]). Also, we shall present the results of the numerical evaluation of eq. (|7.21|) 
for two values of 7. 

(I) Short-range behaviour: r±_ <^ 1/Qs{t) 

To estimate the integral in eq. ( [7.21| ), we divide the integration range into three 
domains separated by the external scales l/r_|_ and Qs{t)- Thus, we need to consider the 
following domains: (I-a) <^ k±, (I-b) (5s(t) ^ ^ l/r_|_, and (I-c) k± <^ <5s(t). 

In domain (I-a), kj_ is the largest scale in the problem, so kj_-r± ^ 1 and Q'i{T)/k'j_ -C 
1. Then, one can ignore the rapidly oscillating exponential e^^-L-^±^ and also expand the 



31 



0.01 0.1 1 10 100 1000 

ki / Qs(xo) 



Figure 5: Energy dependence of <^r(^±) for 7 = 0.64. We have plotted V5r(^±) as a 
function of k^/Qsijo) (with tq some value of reference) for six values of r. The lines, 
from the bottom to the top, correspond successively to r = tq, tq + 2, ■ ■ ■ , tq + 10. The 
increase with r is exponential at high momenta (giving equidistant curves in this log-log 
plot), but only logarithmic at low momenta. 
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logarithm as ln(l + x) ~ x. This gives: 

In domain (I-b), one can expand both the logarithm, as above, and the exponential, 
and keep only the leading term which survives after the angular integral. That is, we can 
replace 1 — e^^-^'^^ i^l^l; thus obtain: 

with a = l/[4c7(l — 7)]. Note that the case 7 = 1 requires a separate treatment. 

Lastly, in domain (I-c), we can expand the exponential as above, and approximate 
ln(l + x) ~ Inx. We obtain: 

fi(i-) ^ 1 ri dkl In = 1 rlQlir). (7.24) 

Putting the results altogether, it is clear that, when 7 = 1, the leading contribution 
comes from domain (I-b), since this is enhanced by a large logarithm. This gives: 

exp{-l.MMl„-jJ^}, (7.25) 

which is formally similar to the MV model (compare to eq. ([4.8|) ), with the crucial dif- 
ference, however, that the scale in the logarithm is now set by the (hard) saturation 
momentum Qs(r), rather than by Aqcd- This reflects the fact that, due to colour neu- 
trality at saturation, the quantum effective theory is free of infrared problems. 

Note also that the exponent in eq. ( [7.25| ) vanishes at r±_ = 1/Qs{t), which indicates 



that this expression is valid only at very short distances, r± -C 1/Qs{t), and thus cannot 
be used to study the approach towards saturation. For this latter purpose, a more refined 
calculation is necessary, and for 7 = 1, this can be even performed analytically. We shall 
return to this calculation towards the end of the section. 

When 7 < 1, on the other hand, it is clear from eqs. (|7.22|) and ( |7.23|) that domains 
(I-a) and (I-b) contribute on the same footing, and give together the leading behaviour. 
This shows that the integral in eq. ( |7.21| ) is now dominated by momenta k± ~ l/''"±, as 



expected for the BFKL dynamics. We thus deduce that, for r± -C 1/Qs{t) (but within 
the scahng window ( p.ll|) ), the scattering amplitude has the right scaling form: 



^ ^ {rlQl{r)y (7.26) 

7<1 



(with K some unknown number) in agreement with eq. ( [6.121 ). 



Note that k± is the momentum transferred between the incoming dipole, of transverse 
size r_L, and the colour sources that it scatters off in the hadronic target. Thus, the 
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previous discussion shows that a small dipole {r± <^ 1/Qs{t)) undergoes mostly hard 
scattering {k± ^ Qs{t)), that is, it interacts predominantly with the colour sources in its 
vicinity in the impact parameter space. This agrees with the conclusions in Ref. |p. 

(II) Long distance behaviour: r± ^ 1/Qs(t) 

The three domains that we have to consider now are: (Il-a) kj_ S> Qs{t), (H^b) 
Qs{'t) S> k± 3> l/r±, and (II-c) l/r± ^ Their respective contributions are easily 
estimated, with the following results: 



Q(II-b) _ i /"^'^^^ ( Q'ii^) \ _ }_( 1 2^2/ ^ 



(7.28) 

^ ijll ^ dkl\n ^±-\nrlQl{r). (7.29) 

Clearly, to leading-log accuracy, the dominant contribution comes from domain (Il-b), 
and gives: 

^.(r^) ocexp|-l(lnriQKr))'| , (7.30) 

in agreement with the results in Refs. [^, |^, |10|, |Tl|]. In this case, the dipole scattering is 
dominated by relatively long-range interactions (fc_L ^ Qs{t)i with k^ ^ 1/?"^ though), 
that is, by the interactions with the saturated gluons. 

Eq. (|7.30|) describes the approach of the S'-matrix element towards the "black disk" 
limit Sr{r±) = in the regime where the dipole is already very large. Thus, when this 
functional form starts to apply, ST-{r_C) is already significantly smaller than one (i.e., the 
proportionality coefficient in eq. (|7.3CI| ) should be rather small). 

Let us finally return to the case 7 = 1, and show that, in this case, the behaviour of 
Sr{r±) around the saturation scale can be studied analytically. To this aim, it is preferable 
to return to the expression of Sr{ri_) in eq. ( |3.19|) [with Xy{k_[_) as given by eq. ( |6.13|) ], 
and perform first the momentum integration there, for fixed y. This generates a modified 
Bessel function (see the Appendix for more details), which is then expanded in an infinite 
series, and the result is integrated over y term by term. The final result reads: 

By using this exact series expansion, one can study the approach towards saturation from 
shorter distances. Namely, for r_|_ -C 1/Qs{t), the dominant contribution comes from the 
n = 1 term in eq. (|7.31|) , which yields 

1 



7=1 I 4c 



ln^^+(l + 2^(2) + 21n2)j}. (7.32) 
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1.2 



Interpolation with y = 0.64 
Short distance approximation 
Long distance approximation 



0.8 



0.6 



0.4 




r±Qs(x) 

Figure 6: The scattering amplitude 1 — Sr{r±) for 7 = 0.64 plotted as a function of 



fi.Qsi.T')- Red solid line: numerical evaluation of eq. (|7.21| ); green dashed line: short dis- 
tance approximation, eq. (|7.26| ); blue dotted line: long distance approximation, eq. ( |7.30|) . 



As compared to the leading-log result in eq. ( [7.25|) , the expression above displays also the 
subleading term under the logarithm. Clearly, when r±_ approaches 1/Qs{t) from below, 
this formally "subleading" term gives the dominant behaviour. 

In Figs. ^ and [7|, we represent the results of the numerical evaluation of the scattering 
amplitude Mriri) = 1 — >S'T-(r_L), eq. (|7.21|) , for two values of the anomalous dimension: 



7 = 0.64 (Fig. ^) and 7 = 1 (Fig. 0). In the same plots, we display, for comparison, 
the corresponding predictions of the various approximations derived previously in this 
subsection. 



8 Conclusions and perspectives 

In this paper, we have proposed a simple Gaussian approximation for the weight function 
for the effective Colour Glass description of the small-x gluons in the hadron wavefunction. 
The main virtue of this approximation is its simplicity. The new weight function is as 
simple to use in practice as that of the original McLerran-Venugopalan model, but it 
improves over the latter by including the correlations among the colour sources associated 
with quantum evolution towards small x and gluon saturation. This ensures, in particular, 
the removal of the infrared divergences inherent to the MV model, without invoking non- 
perturbative physics: Gauge-invariant observables come out infrared finite because of the 
colour neutrality of the saturated gluons over the (relatively short) scale 1/Qs{t). More 
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1.2 



Interpolation with y = 1 
Short dist. approx. (leading) 
Short dist. approx. (next leading) 
Long distance approximation 




riQs(x) 



Figure 7: The scattering amplitude 1 
Red solid line: numerical evaluation of eq 



Sr{r±) for 7 



1 plotted function of r±Qs{T). 
green dashed line: leading-order short 



distance approximation, eq. (|7.25| ); blue dotted line: next-to-leading-order short distance 
approximation, eq. ( |7.32| ); magenta dotted line: long distance approximation, eq. ( |7.30| ). 



generally, this Gaussian weight function implements the BFKL dynamics (for both the 
gluon distribution, and the dipole-hadron scattering) at high transverse momenta, and 
gluon saturation at low transverse momenta, with the change between these two regimes 
taking place at the saturation scale, as it should. 

In its simplest form, the kernel of the Gaussian is a scaling function, eq. ( |6.13D , which 
is strictly valid up to momenta k± ^ Q^r)/ Aqcd, and involves three dimensionless 
parameters: the anomalous dimension 7, the logarithmic derivative of the saturation mo- 
mentum A = c? In Q^(r)/(ir (throughout this paper, we have rather denoted this quantity 
as A = cds ; see, e.g., eq. (|6.7D), and the initial rapidity tq = ln(l/xo) which appears when 



the saturation momentum, eq. (|6.7|), is rewritten as Ql{x) = Q^e 



)2pA(T-ro) 



Ql{xo/x)^, with 



Qo some convenient scale of reference (e.g., Qo = 1 GeV). Of course, within the present 
approximations, the values of 7 and A are determined by the leading-order BFKL dy- 
namics. But since we expect geometric scaling and the exponential rise of the saturation 
scale with r to be more general than just LO BFKL, we can also treat these quantities as 
free parameters, to be fitted from the data, or taken from some more refined calculations. 
For instance, the fits to the electron-proton deep inelastic scattering data at HERA using 
the "saturation model" [^, ^ provide a value A ~ 0.3, which, remarkably, appears to 



be consistent with a recent calculation |2C] of the saturation scale using the NLO BFKL 
dynamics (with the RG-improved kernel proposed in Refs. fS^). 



The Gaussian weight function that we have constructed and justified here can be 
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immediately used to supplant the MV model in its applications to the physics at small x. 

For the problems involving only one high-density hadronic system (one large nucleu^, 
or a very energetic hadron), the calculations can be performed quasi-analytically, since 
the solution to the corresponding classical field equations is known in analytic form (cf. 



Section 2.1). This includes various studies of proton-nucleus (pA) collisions |§, ^ 
or of ultra-peripheral nucleus-nucleus {AA) collisions |5^. Also, the expression for the 
dipole-hadron scattering amplitude that we have obtained in this context, eq. ( [7.21|) , can 
be used in analyses of the DIS data at HERA, as an alternative to the more phenomeno- 
logical parametrizations proposed in Refs. ^9|, 0. It would be interesting, in this respect, 
to compare eq. ( |7.21| ) with the numerical solution of the BK equation p2| , |23| , or 
with the results of the lattice simulations of the full RGE [5^] (the latter should be soon 
available We expect our analytic formula ( |7.21| ) to be numerically close to such 

more complete calculations, while being much simpler to use in practice, and allowing 
for the freedom to choose (or fit) the values of the parameters 7 and A (in the numerical 
calculations alluded to above, these parameters are rather fixed to their respective BFKL 
values, by the dynamics of the relevant equations). 

For applications of the CGC effective theory to heavy ion collisions, the classical 
field equations can be numerically solved on a lattice, as in Refs. p2i, and their results 



averaged over the colour charge distributions in the incoming nuclei, with the Gaussian 
weight function. As compared to the MV model, the use of the new Gaussian in such 
applications would eliminate the sensitivity to the poorly known non-perturbative physics. 
Still within the context of AA collisions, we note that the simple analytic form we have 
derived for the unintegrated gluon distribution, eq. ( |7.1^ ), can be directly implemented 
in analytic calculations of the multiparticle production, like those in Refs. EBl pTl. 
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A Appendix 



In this Appendix, we derive the analytic expression ( [7.31j ) for the exponent of the S- 
matrix element, which refers to eq. ( [7.21| ) with 7 = 1. We first perform the momentum 
integration in eq. ( p.l9| ), where Xy{k±) is now replaced by its interpolation in eq. ( |6.13|) 
(with 7 = 1). It is useful to decompose this latter as follows: 



7=1 



1 



1 



(A.l) 



When this function is inserted in the integrand of eq. ( p.l9| ), its product with 1 — e^'^^'^-L is 
naturally decomposed into four terms. Since the overall integral is finite, we are allowed 
to evaluate these four terms separately, provided we use the same regularization scheme. 
All these terms are represented by a simple integral (after angular integration, r = |r_|_|): 



dk Jo{kr) 



k 



KoiPr) 



J.2 + p2 

where Ko{z) is the modified Bessel function (with ip^z) = dliaT{z)/dz) 



(A.2) 



K,(z) = - 



lnz + (ln2-7)+$^-^(-) (^(n+l)-ln 

n=l ^ ' 



(A.3) 



For example, the term J{d'^k^/k\), which has both ultraviolet and infrared divergences, 
can be regularized as / 2T[kdk jQ{k/ K)/{k'^ + yU^), with A and /i being the UV and IR 
cutoffs, respectively, and thus brought into the form of eq. (|A.2| ). When adding the 
contributions of the four terms, all divergences cancel out, as they should, and, in the 
limit A — i> 00 and /i — >^ 0, we are left with the following result: 



^ = 2a. j dy | ln(a(y)r) - (ln2 - 7) + Ko{Qs{j)r)]. 



(A.4) 



The integral over the rapidity y can be now performed, by using the energy dependence 
of Qs(j) given in eq. (§]^ : 



7=1 



4 

c 



I ^{lnC-(ln2-7)+i^o(C)} [ 



C = Qs{j)r ] 



C 



(A.5) 



In the last line, we have used eq. ( [A.3| ). By integrating this series term by term, one 
obtains the final result in eq. ( [7.31| ). 
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